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We extend the analysis of [Zhou and Sornette, Physica D 165, 94-125, 2002] showing statistically 
significant log-periodic corrections to scaling in the moments of the energy dissipation rate in experi- 
ments at high Reynolds number (~ 2500) of three-dimensional fully developed turbulence. First, we 
develop a simple variant of the canonical averaging method using a rephasing scheme between differ- 
ent samples based on pairwise correlations that confirms Zhou and Sornette's previous results. The 
second analysis uses a simpler local spectral approach and then performs averages over many local 
spectra. This yields stronger evidence of the existence of underlying log-periodic undulations, with 
the detection of more than 20 harmonics of a fundamental logarithmic frequency / = 1.434 ± 0.007 
corresponding to the preferred scaling ratio 7 = 2.008 ± 0.006. 
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I . I. INTRODUCTION 

The small-scale physical quantity that has received the most attention in the context of a phenomenological de- 
scription of hydrodynamical turbulence is the dissipation rate of kinetic energy ffl. Richardson's picture of turbulent 
cascades in which large eddies break down into smaller ones receiving a certain fraction of the flux of kinetic 
energy from larger scales, is a multiplicative process. The hypothesis that eddies of any generation are space filling 
. led to the famous K41 theory B 0|. A conceptually appealing view, dating back to Obukhov Q and Kolmogorov 
' ||, visualizes the transfer of kinetic energy as a self-similar multiplicative process. This view is still the point of 
departure of many phenomenological intermittency models today. The thought of the general theory of multifractals 
was proposed by Mandelbrot in the study of certain high-order moments with a general cascade model of random 
curling. A phenomenological model of intermittency, referred to as the /3-model ||, analyzed the multiplicative cas- 
cade of the breakup of eddies without space filling leading to an energy transfer to a smaller and smaller fraction of 
the total space. Benzi, Paladin, Parisi & Vulpiani |J and Frisch & Parisi ]l0) have argued for a singular structure 
and for a multifractal nature of the energy dissipation in three-dimensional fully-developed turbulence. Meneveau & 
Sreenivasan fi"i"|| presented a simple multifractal model, the binomial model for fully developed turbulence based on 
the two-scale Cantor set, which is known as the p-model. This model and the specific multifractal facets of turbulent 
energy dissipation were investigated experimentally by [[l2|, using hot wire measurement technique. These works 
clarified the phenomenological nature of the multifractality of turbulent energy dissipation. In particular, the view of 
the energy cascade as belonging to the general class of multiplicative processes opens up the possibility that the scale 
invariant properties reported for turbulence flows might be broken in part into a discrete scale invariance, leading to 
the observation of log-periodic corrections to scaling |l3| . 

Based on theoretical argument and experimental evidence, Sornette |l4j conjectured that structure functions of 
turbulent time series may exhibit log-periodic modulations decorating their power law dependence. He stressed the 
need for novel methods of averaging and proposed a novel "canonical" averaging scheme for the analysis of structure 
factors of turbulent flows in order to provide convincing experimental evidence. The strategy was proposed to 
determine the sample-specific control parameter r c and translate it into a specific "phase" in the logarithm of the 
scale lnr which, when used as the origin, allows one to rephase the different measurements of a structure factor in 
different realizations. 
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A first systematic search of discrete scale invariance in turbulence was carried out in freely decaying two-dimensional 
turbulence |T^| . The number of vortices, their radius and separation were found to display log-periodic oscillations as a 
function of time with an average log-frequency of ~ 4 — 5 corresponding to a preferred scaling ratio of ~ 1.2 — 1.3. Most 
recently, significant evidence for the presence of intermittent discretely self-similar cascades was obtained from the 
analysis of the moments of the energy dissipation rate in three-dimensional fully developed turbulence JT(| . The results 
of spectral analyses of the canonically averaged local scaling exponents of high moments of the energy dissipation rate 
obtained in p6[ are summarized in Fig. [lj The fundamental log- frequency is estimated to be / « 1.44 which suggests 
a preferred scaling ratio of 7 = e 1 ^ rj 2. In a simple multifractal model was introduced to explain the fact that 
the log-frequency is independent of the order of the moments. 

As was explained in , the difficulty in detecting log-periodic oscillations (the hallmark of discrete scale invariance) 
in various functions such as moments of the energy dissipation rate stems from the fact that, while the physics of the 
cascade may be universal with an universal scaling factor 7, the nucleating scales (translating in phases of reference 
in the log-periodicity analysis) have to change from one realization to another as the cascade may be dynamically 
triggered from a variety of scales £ a from the integral scale and deep in the inertial range. Therefore, the corresponding 
phases are expected to be non-universal. Averaging a signal over several such transient cascades will thus wash out 
the information on log-periodicity by "destructive interferences" of the log-periodic oscillations, giving them the 
appearance of noise. 

The previous analyses of real 2-D and 3-D turbulence data addressed this problem by using two versions of the 
canonical averaging approach, devised precisely to rephase different realizations and thus alleviate the problem of the 
random phases. Using n realizations of the data, the first version consists in performing a canonical averaging on 
the n samples to get an average realization and then analyzes its spectral properties; the second version performs a 
spectral analysis of each of the n samples and then averages the n periodograms. In order to perform these averages, 
we need to have different realizations. For the turbulence data analyzed in |l6| and revisited here, a long time series 
of velocity increments is divided in many samples, each lasting about 1/5 of a turn-over time, so that we can avoid 
scrambling of the phases by a suitable canonical averaging as explained below. 

The goal of the present note is to revisit these two schemes and propose two novel and simpler versions. By 
simplifying the analysis and by providing independent implementations, we thus confirm further the presence of log- 
periodicity in the moments of the energy dissipation rate, and thus the existence of discrete scale invariance in 3D 
fully-developed hydrodynamical turbulence. 



II. REPHASING OF LOG-PERIODICITY BY THE CROSS-CORRELATION METHOD 

In the first scheme of canonical averaging, the rephasing operation is the central part. In |lq ], the so-called central 
maxima approach was applied to the detrended logarithm of the moment M q (l/rj) of the dissipation rate as a function 
of the logarithm \n.(l/rj) of the scale £ (in units of the Kolmogorov microscale rj = 0.195 mm [jf6|): it consists in 
finding the maximum of the detrended logarithm of the moment closest to the middle point of a single sample and set 
the logarithm of the corresponding scale (its abscissa) as the origin of the logarithm of the scales. The new abscissa, 
called A, becomes sample dependent and is used to superimpose and then average the detrended moments over the 
different samples. This approach was found to work satisfactorily except that the canonically averaged series present 
decreasing amplitudes with increasing |A|. 

The central maxima rephasing approach used in fl(| determines an absolute sample-dependent phase which is 
specific to each single sample and then rephases all samples. The alternative idea explored here is based on the 
remark that the detection of phases can be performed in relative value, that is, by estimating the phase shift between 
all possible pairs of samples. The cross-correlation function is used to calculate the phase-shift between two samples 
jl7j . This relative phase-shift approach thus uses the calculation of all phase-shifts between pairs of samples to rephase 
all samples by choosing arbitrarily one of them as the reference. 

As in |L6|, we use the high- Reynolds turbulence longitudinal velocity data collected at the SI ONERA wind tunnel 
by the Grenoble group from LEGI jl8). We use a set of 20 records, each containing 320 samples in the inertia range, 
each sample having 2 11 points. These 2 11 points correspond approximately to one-fifth of a turn-over time |l6|. We 
calculate the local moment exponent r (q, \n(£/rj)) as the logarithmic derivative of M q (£/rj) with respect to £/rj: 

r( q ,W/v))= dHe/v) ■ (1) 

We shall investigate in the sequel £/rj G [230,1843] which is well within the inertial range. To obtain r (g, ln(£/j])) 
from (g), we used a smoothing low-pass digital filter, namely, the Savitzky-Golay filter. This filter approximates 
the underlying function in a moving window by means of polynomials of a low degree. The filter coefficients arc 
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choosen so that, unlike other filters, several low-order moments of the filtered function are preserved, which is a very 
desirable property (see for details). The Savitzky-Golay filters are applied to ]n[M q (£/rj)] expressed as a function 



of hx{£/r]) and the local derivative is given from the the analytic derivative of the fitted polynomial |19 . There are 
two parameters in the Savitzky-Golay filters, i.e., the width of the sliding window Nl + Nr + 1 where Nl and Nr are 
the numbers of points to the left and to the right of the investigating point and the order M of fitting polynomials. 
We used symmetric windows such that Nl — Nr. A total of 24 filters were applied to test the robustness of the 
results with respect to the choice of the filters: Nl — 5, • • • ,10 and M = 4, • • • ,7. 

We recall that, in the central maxima approach of ||lq| , one identifies the sample-dependent reduced scale £ c hl 
closest to the central point of the scale-interval at which r (q, ln(£/rj)) is maximum for each sample. In contrast, the 
correlation approach takes an arbitrary sample as the reference sample and evaluates the phase shift measured by the 
cross-correlation function of all other samples relative to this reference. We then rephase all functions r (q,ln(£/rj)) 
of the variable \n(£/rj) so that the cross-correlation function of all sample pairs becomes peaked at zero lag, after the 
rephasing has been performed. 

We then average these rephased functions t (q, A) over the 320 samples. This defines the local average exponent 

D(A) = (r(q,A)) , (2) 

where the average is performed over the different samples at fixed A = ln(£/ry) — hi(£ c /r]) (where \n(£ c /r)) vary from 
sample to sample and \n(£/r]) is adjusted accordingly to ensure the re-phasing). Fig. |] shows four examples of D(A) 
as a function of A based on the correlation rephasing approach for four different filters with Nl — 6 and M = 4, 5, 6, 7 
resp ectively. Notice that the amplitude of the oscillations decreases with decreasing A. This is in line with the concept 
|L6| that, the smaller is the scale £, the more complex is the superposition of intermittent cascades as the number 
of starting scales for the cascade proliferate when the nucleating scale departs from the integral scale. This should 
lead to a stronger self-averaging property, as we observe here when A decreases. Furthermore, since the samples have 
finite sizes, -D(A) is well-defined only for not too large |A| because there are fewer samples involved in the averaging 
for larger |A|. In practice, we will restrict our analysis to IA| smaller than a threshold A m taken here approximately 
equal to 1.0, as suggested from a visual inspection of Fig. 0. 

The next step in the analysis consists in performing a spectral analysis of D(A). We used a kind of spectral analysis 
known under the name of the Lomb spectrum (see pi| for details). It is designed to analyze a non-uniformly sampled 
functions, fitting them to a sine function of varying frequency on a non- uniform grid by means of the least square 
method. This provides the analog of a standard periodogram. We thus obtain 24 such Lomb periodograms for each 
record (of 320 samples), one for each filter. We then repeat the above steps for all the other 19 records of of 320 
samples each and we thus obtain a total of 480 Lomb periodograms. Alternatively, we have grouped the 20 records 
into a single large time series with 6400 samples on which the canonical averaging is performed directly. This increased 
statistics improves the significance level of the Lomb peaks in the resultant periodograms. The comparison of our 
method on different partitions of the data allows us to test the stability and robustness of the results. 

Having computed the 480 Lomb periodograms (corresponding to 24 Lomb periodograms for each of the 20 records), 
we perform an average over all these 480 periodograms. This analysis was carried out for three different choices of the 
threshold A m = 1.0,0.8,0.6. The averaged Lomb powers are shown in Fig. ||. All the three curves have significant 
peaks at f% — 2/ln2 and — 3/ln2, which suggest a fundamental log-frequency / = l/ln2. The periodogram 
with A m = 1.0 also possesses many spurious weaker peaks. The periodogram with A m = 0.8 exhibits the highest 
Lomb peaks and clearest harmonics. The difference in peak heights for A m — 0.8 compared with A m — 0.6 is 
partially explained by the fact that the former has more data points. It is indeed known that the Lomb power is 
approximately proportional to the number of analyzed points |20| . The results summarized in Figure [| thus provide 
a strong confirmation of the analysis of pM on the presence of discrete scale invariance. 



III. DIRECT DETECTION OF LOG-PERIODICITY BY FAST FOURIER TRANSFORM 

We now turn to the second method of implementation of the canonical averaging discussed in the introduction. In 
a nutshell, the methods consists in performing directly a spectral analysis of each sample and then average over the 
obtained spectra. In order to improve the signal-to-noise, we first filter the high-frequencies of the function \n[M q (£/rj)] 
of the variable \n(£/rj) by performing an interpolation of the slightly unevenly sampled third-order moment values 
M3 (£/??) with cubic splines. This provides us with equally sampled points that are required for the use of a fast Fourier 
transform (FFT). We then define the first difference and subtract the mean which defines the function 8\n[M^(£/'q)] 
The periodogram of <51n[M3 (£/??)] for each sample is determined via a fast Fourier transform. We then perform an 
averaging over all the periodograms. 

Specifically, the calculation of the third-order moment M^{£/rf) is done as in Sec. II with the same data sets. We 
use one record with 320 samples. The number of the evenly sampled points generated by the spline interpolation is 
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257, such that we will have 256 points at which the function 5\n[M 3 {£/vi)\ is sampled. Fig. || shows a realization of 
M 3 (£) for one specific sample. For this, 6i\n[M 3 (£/r])] is obtained as 

5i \n[M 3 {i/r,)] = ln[M 3 (W»?)] ~ HM 3 (ii/v)] ~ M, (3) 

where i = 1, • • ■ ,256 and /i = (ln[M3(£j+±/rf)] — \n[M^{£j/r})\)j, The resulting detrended third-order moment is 
shown in Fig. [s] and its periodogram is depicted in Fig. The vertical lines in Fig. |^ indicate log-frequencies 
equal to fc/ In 2 with k — 1, 2, ■ • ■ ,21. One can observe well-defined peaks close to the frequencies corresponding to 
k = 1,2,3,4,5,6,9,12,15,16,17,21. However, there are also peaks in between these vertical lines. Therefore, the 
signature of log-periodicity in Figs. |^-|| is ambiguous. 

We then average over the 320 periodograms, one for each sample, and get Fig. ^. Most peaks at log-frequencies other 
than kj In 2 have been washed out by the averaging procedure. The vertical dashed lines correspond to log-frequencies 
fk = k/ In 2 with k — 1, • • • ,21, which are the harmonics of a fundamental log- frequency l/ln2 corresponding to a 
preferred scaling ratio equal to 7 = 2. The local peaks defined preferred frequencies F m that are shown as open 
circles. Most of the peaks (m = 3, • ■ • ,18, 20, 22, 24, 27) in Fig. are very close to the vertical dotted lines which can 
be interpreted as harmonics of / = 1/ In 2. The background noise has a bell-shape. Comparing Fig. |?j with Fig. |6[ the 
averaging of the spectra over the 320 samples has suppressed peaks that can be interpreted as noise and has improved 
remarkably the significance level of the other Lomb peaks. 

In order to estimate as accurately as possible the numerical value of the fundamental log-frequency, Fig. || plots with 
open circles the log-frequencies F m of all the peaks found in Fig. [7] as a function of the peak sequence number m. One 
observes a straight line with slope / = 1.420 covering 16 harmonics. This value is very close to the theoretic ansatz 
/ = l/ln2 w 1.4427. An estimation of the uncertainty in the determination of / is obtained by removing one point 
in each data set, and redoing the whole procedure. Then, the standard deviation calculated over these samples with 
one removed point provides a quantification of the error in the determination of /. In order to include more points in 
the fit, we collected peaks near the vertical lines in Fig. [?] with log-frequencies fk and plotted fk with open squares 
in Fig. H as well. Note that fk with fc = 1, • • • , 20 corresponds to peaks enumerated m = 3, • • • , 18, 20, 22, 24, 27 
mentioned above. Again, a very good linear fit is obtained with slope / = 1.427, which provides a slightly better 
estimate of 1/ In 7. We estimated the slopes of 19 other experimental records which, together with the slope from the 
first record, give the estimate of the fundamental log- frequency / = 1.434 ± 0.007. An alternative estimate of / is 
the mean of the ratios fk/k. This gives / = 1.455 ± 0.023. In this estimation, the first point (fc = 1, //. = 1.864) was 
discarded. 

Overall, we conclude that this new analysis agrees with and strengthens strongly previous results obtained in Sec. II 
and in |l6| on the observation of log-periodicity. It is also quite surprising and somewhat gratifying to observe so 
many harmonics. 



IV. CONCLUSION 



In this paper, we have provided two novel analyses for the search of log-periodic corrections to scaling in the moments 
of energy dissipation of three-dimensional fully developed turbulence. The first analysis introduces a pairwise cross- 
correlation approach for the determination of a phase-shift, which is an alternative of the central maxima rephasing 
approach developed in [ fl6| and is more natural and straightforward to implement. The resulting averaged Lomb 
periodogram have two significant peaks near log- frequencies equal to 2/ In 2 and 3/ In 2 and low peaks near other 
fc/ln2 with k = 1,4, 5. We have checked that these results are robust when varying the range of scales over which 
the analysis is performed. 

The second method uses a more straightforward sample-by-sample spectral analysis, with a simple cubic spline 
interpolation (rather than the more delicate Savitzky-Golay filters) and a standard FFT. Averaging over the hundreds 
of spectra provides a strong log-periodic signals with more than 20 harmonics of what appears to be a fundamental 
logarithmic frequency / = 1.434 ± 0.007. The corresponding preferred scaling ratio is thus estimated to be 7 = 
2.008 ± 0.006. This novel analysis lifts possible remaining worries that one could keep that the results of |16 could 
be somehow due to the use of the Savitzky-Golay filters. 

The present work further confirms the presence of log-periodicity in the energy transfer of fully developed turbulence, 
as well as provide additional evidence of the presence of multifractality in the energy dissipation p2[| . However, the 
preferred scaling ratio 7 sa 2 found here does not necessarily imply that the correct model for turbulence is the 
binomial multifractal cascade plf . We ascribe our results as an imprint of intermittent cascades from discrete 
hierarchical dissipation in turbulence |lrJ . 
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FIG. 1: Average of 2400 (continuous line) (resp. 480 (dashed line)) Lomb periodograms corresponding to 24 Savitzky-Golay 
niters per sample. The vertical dashed lines correspond to log-frequencies equal respectively to fk = fc/ln2 with k = 1, ••• ,5. 
These values correspond to increasing harmonics of a putative fundamental frequency /i = 1/ In 7 associated with the scaling 
ratio 7 = 2. (See Ref. B for details.) 
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FIG. 2: Four examples of the canonically averaged variable D(A) as a function of the rephased log-scale A based on the 
correlation rephasing approach. The amplitude of the oscillations increases with increasing A. See text for explanations. 
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FIG. 3: Average over all 480 Lomb periodograms corresponding to 24 niters per sample. The vertical dashed lines correspond 
to log-frequencies equal to fk = k/ In 2 with k — 1, • ■ • ,5. These values correspond to increasing harmonics of a fundamental 
frequency /i = 1/ In 7 associated with the scaling ratio close to 7 = 2. 
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FIG. 5: The zero-mean detrended <51n(M3) as a function of ln(£/rj) for the spline interpolation shown in Fig. 
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FIG. 6: Fourier periodogram of S ln(Afs) shown in Fig. [s| The vertical lines indicate log-frequencies k/ In 2 with k = 1, • • • , 21. 
There are peaks close to the vertical lines with k = 1,2, 3, 4, 5, 6, 9, 12, 15, ... and also peaks that can not be integrated as integer 
multiples of l/ln2. 
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FIG. 7: Average of all Fourier periodograms over 320 samples. The vertical dashed lines correspond to log-frequencies equal 
to fk = fc/ln2 with k = 1, • • • , 21. The peaks are indicated with open circles. Most of the peaks are very close to the vertical 
dotted lines which can be interpreted as harmonics of / = l/ln2. 




FIG. 8: Determination of the fundamental log-frequency with a linear fit: the open circles are all peaks indicated in Fig. [7j 
and the open squares are those close to the vertical lines in Fig. (?|. Excellent linear fits are obtained with the slopes 1.420 (o) 
and 1.427 (□), close to the value l/m2 = 1.4427. 



